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Abstract. Let G be a bipartite graph with positive integer weights on the edges and without 
isolated nodes. Let n, and W be the node count, the largest edge weight and the total weight of 
G. Let k{x,y) be log x/ log(a::-^/y). We present a new decomposition theorem for maximum weight 
bipartite matchings and use it to design an 0{y/nW/k{n,W/N))-time algorithm for computing a 
maximum weight matching of G. This algorithm bridges a long-standing gap between the best 
known time complexity of computing a maximum weight matching and that of computing a maximum 
cardinality matching. Given G and a maximum weight matching of G, we can further compute the 
weight of a maximum weight matching of G — {u} for all nodes u in 0{W) time. 
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1. Introduction. Let G ~ {X,Y,E) be a bipartite graph with positive inte- 
ger weights on the edges. A matching of G is a subset of node-disjoint edges of G. 
Let mwm(G) (respectively, mm(G)) denote the maximum weight (respectively, car- 
dinality) of any matching of G. A maximum weight matching is one whose weight is 
mwm(G). Let N be the largest weight of any edge. Let W be the total weight of G. 
Let n and m be the numbers of nodes and edges of G; to avoid triviality, we maintain 
m = f2(n) throughout the paper. 

The problem of finding a maximum weight matching of a given G has a rich his- 
tory. The first known polynomial-time algorithm is the 0(n'^ )-time Hungarian method 
[p^ . Fredman and Tarjan used Fibonacci heaps to improve the time to 0{n{m + 
nlogn)). Gabow Q introduced scaling to solve the problem in 0(n'^/^m log A^) time 
by taking advantage of the integrality of edge weights. Gabow and Tarjan |7| im- 
proved the scaling method to further reduce the time to 0{^/nm,\og{nN)). For the 
case where the edges all have weight 1, i.e., N = \ (and W = m), Hopcroft and Karp 
[pl gave an 0{^/TlW)-time algorithm, and Feder and Motwani 0| improved the time 
complexity to 0{^/nW/k{n,m,)), where k{x,y) = \ogx/ \og{x^ /y). It has remained 
open whether the gap between the running times of the Gabow- Tarjan algorithm and 
the latter two algorithms can be closed for the case where W = o(m log(nA^)). 

We resolve this open problem in the afhrmative by giving an 0{y/nW/k{n, W/N))- 
time algorithm for general W . Note that W/N — m when all the edges have the same 
weight. The algorithm does not use scaling but instead employs a novel decomposition 



theorem for weighted bipartite matchings (Theorem 2.2). We also use the theorem to 
solve the all-cavity maximum weight matching problem which, given G and a maxi- 
mum weight matching of G, asks for mwm(G — {m}) for all nodes u in G. This problem 
has applications to tree comparisons 1^,0. The case where = 1 has been studied 
by Chung Recently, Kao, Lam, Sung, and Ting ||l2| gave an 0{^Jnm\ogN)-\,ime 
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(a)G (b)G, (c)C, (d)G^ 

Fig. 2.1. Consider h = 1. G is decomposed into Gii and G^; Ch is a minimum weight cover 
ofGh- 



algorithm for general N . This paper presents a new algorithm that runs in 0{W) 
time. 

Section ^ presents the decomposition theorem and uses it to compute the weight 
of a maximum weight matching. Section |^ gives an algorithm to construct a maximum 
weight matching. Section ^ solves the all-cavity matching problem. 



2. The decomposition theorem. In |2.1| , we state the decomposition theo- 
rem and use the theorem to des ign an algorithm to compute the weight mwm(G) in 
0{'s/nW/k{n,W/N)) time. In §2^, we prove the decomposition theorem. In we 
further construct a maximum weight matching itself within the same time bound. 

2.1. An algorithm for computing mwm(G). Let V{G) be the node set of 
G, i.e., X U y. Let w{u^ v) denote the weight of an edge uv G G; if u is not adjacent 
to V, let w{u,v) = 0. A cover of G is a function G : X UY {0, 1,2, . . .} such 
that C{x) + C{y) > w{x,y) for all x € A and y &Y. Let w{C) — J2zexuY ^(•^) 
the weight of G. G is a minimum weight cover if w{C) is the smallest possible. Let 
mwc(G) denote the weight of a minimum weight cover of G. A minimum weight cover 
is a dual of a maximum weight matching as stated in the next fact. 

Fact 2.1 (see ]l|). Let C he a cover and M be a matching of G. The following 
statements are equivalent. 

1. C is a minimum weight cover and M is a maximum weight matching of G. 

^- T,uv&mW{u,v) =Y,u&XUYCi'^)■ 
3. Every node in {u \ G{u) > 0} is matched by some edge in M, and C{u) + 

C(v) = w(u,v) for all uv S M. 

For an integer h £ [1, A], we divide G into two hghter bipartite graphs Gh and 

G^ as follows: 

• Gh is formed by the edges uv of G with w{u, v) S [N — h + 1, N]. Each edge 
uv in Gh has weight w{u,v) — (A — h). For example, Gi is formed by the 
heaviest edges of G, and the weight of each edge is exactly one. 

• Let Ch be a minimum weight cover of Gh ■ G^ is formed by the edges uv of G 
with w{u, v) — Ch{u) — Gh {v) > 0. The weight of uv is w{u, v) — Ch{u) — Ch{v)- 

An example is depicted in Figure 2.1. Note that the total weight of Gh and G^ is at 
most W. 

The next theorem is the decomposition theorem. 

Theorem 2.2. mwm(G) = mwm(G;i) + mwm(G^); in particular, mwm(G) = 
mm(Gi) -I- mwm(Gf ). 
Proof SeejO. □ 

Theorem 2.2 suggests the following recursive algorithm to compute mwm(G). 
Procedure Compute-MWM(G) 
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1. Construct Gi from G. 

2. Compute mm(Gi) and find a minimum weight cover Ci of Gi. 

3. Construct Gf from G and Gi. 

4. If Gf is empty, then return mm(Gi); otherwise, return mm(Gi)+Compute- 
MWM(Gf). 

Theorem 2.3. Compute-MWM(G) /inds mwm(G) inO{^/nW/k{n,W/N)) time. 



Proof. The correctness of Compute-MWM follows from Theorem 2.2. Below, we 
analyze the running time. We initialize a maximum heap |^ in 0(m) time to store 
the edges of G according to their weights. Let T{n, W, N) be the running time of 
Compute-MWM excluding this initialization. Let L be the set of the heaviest edges 
in G. Then Step 1 takes 0(|L|logm) time. In Step 2, we can compute mm(Gi) in 
0{^\L\/k{n,\L\)) time [|. From this matching, Gi can be found in 0(|-Zj|) time 
[Q. Let Li be the set of the edges of G adjacent to some node u with Gi(m) > 0; 
i.e., Li consists of the edges of G whose weights are reduced in Gf . Let £i = \Li\. 
Step 3 updates every edge of Li in the heap in 0{£i logm) time. As L C Li, Steps 1 
to 3 altogether use 0{y^£i/k{n,li)) time. Since the total weight of Gf is at most 
W — £i, Step 4 uses at most T{n, W — ii,N') time, where N' < N is the maximum 
edge weight of Gf . In summary, for some positive integer £i < W, 

T{n, W, N) = 0{y^£i/k{n, £1)) + T{n, W - t^N'), 

where T{n,0, N') — 0. By recursion, for some positive integers £i,£2, . . . ,£p with 
p< A^andEi<.<p^^ = W^, 



T{n,W,N)^0[ ^ + ^2 



O 



k{n,£i) k{n,£2) 



k{n,£p)^ 



logn 



l<i<p l<i<p ' I 



Since cclogx is convex, by Jensen's Inequality ||lC 



i<i<p 



w 



l<i<p 



Therefore, 



T(n, W,N) = 0\ 



O 



9 w 

Wlogn'' -Wlog- 
iogn V A* 



logn/log(n2/-jv. 



= 0{V^W/k{n,W/N)). 



□ 



2.2. Proof of Theorem 2.2. This section proves the statement that mwm(G) = 
mwm(G;i)-|-mwm(G)^), where G^ is defined according to an arbitrary minimum weight 
cover Gh of Gh. By Fact it suffices to prove mwc(G) — w{Ch) + mwc(G)^). 

To show the direction mwc(G) < w{Gh) + mwc(G^), note that any cover D of 
G^ augmented with Ch gives a cover G of G, where C{u) — Gh{u) + D{u) for each 
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node u of G. Then C{u) + C{v) > w{u,v) for all edges of G. Thus, mwc(G') < 
w{Ch)+nmc{Gt). 

To show the direction w{Gh) + mwc(G^) < mwc(C), let C be a minimum weight 



cover of G. A node u of G is called bad if C(u) < Ch{u). Lemma 2.4 below shows 
that G must have a minimum weight cover C allowing no bad node. Then we can 
construct a cover D of G^ as follows. For each node u of G, define D{u) = C{u) — 
Gh{u), which must be at least 0. D is a cover of G^ because for any edge uv of 
G^, D{u) + D{v) = C{u) + C(v) - Ch{u) - Ch{v) > w{u, v) - G,,(u) - Chiv). Note 
that = w{G) — w{Gh)- Thus, mwc(G^) < w{C) — wiCh), or equivalcntly, 

mwc(G^) + w{Ch) < mwc(G). 

The next lemma concludes the proof of Theorem |2.2| . 

Lemma 2.4. There exists a minimum weight cover of G such that no node of G 
is bad. 

Proof. Suppose, for the sake of contradiction, that every minimum weight cover 
allows some bad node. Then we can obtain a contradiction by constructing another 
minimum weight cover with no bad node. 

Let G be a minimum weight cover of G with u as a bad node, i.e., G{u) < 
Gh(u). Recall that Ch is a minimum weight cover of Gh- Consider a maximum weight 
matching M of Gh- By Fact |]^, since Ch{u) > G(u) > 0, m is matched by an edge in 
M, say, to a node w, and Ch{u) + Gh{v) ~ w{u, v) — {N — h). We call v the mate of u. 
Note that w cannot be a bad node; otherwise, C{u)+C{v) < w{u,v) — {N—h) < w{u,v) 
and a contradiction occurs. 

Since G is a cover of G, C{u) + C{v) > w{u,v). Thus, G{v) > w{u,v) — G(u) > 
N ~ h + Ch(u) + Ch{v) — G{u). Define another cover G' of G as follows. For each 
bad node defined by G, let be the mate of u, define G'(u) = Gh{u) and G'{v) = 
G{v) — {Ch{u) — G(v)). Note that u is not a bad node with respect to G', and neither 
is V since G' {v) > N — h + Gh{v) > Gh{v). For all other nodes x, G'{x) is the same as 
G{x). Therefore, if G' is a cover of G, G' allows no bad node. Also, w{G') = w{G). 

It remains to prove that G' is a cover of G. By the definition of G', G'{v) < G{v) 
if and only if v is the mate of a bad node with respect to G. Suppose G' is not a 
cover of G. Then there exists an edge vt such that G'{v) + G'{t) < w{v,t) and v is 
the mate of a bad node. Recall that the latter implies that G'{v) > N — h + Gh{v). 
In other words, 

G'{t) < w{v, t) - G'{v) < w{v, t)-{N-h)- Gh{v). 

We can derive a contradiction as follows. 

Case 1: w{v,t) < N - h. Then G'{t) < -Gh{v) < 0, which contradicts that 
G'it) > Gh{t) > 0. 

Case 2: w{v,t) > N ~ h. Then Gh contains the edge vt and Gh{v) + Ch{t) > 
wiv,t)-{N -h). Thus, C'{t) < w{v,t) - {N - h) - Gh{v) < Gh{t), which contradicts 
the fact that G' allows no bad node. 

In conclusion, G' is a cover of G. Together with the fact that w{G) = w{C'), we 
obtain the desired contradiction that G' is a minimum weight cover of G with no bad 



node. Lemma 2.4 follows. □ 



3. Construct a mELximum weight matching. The algorithm in §2.1 only 
computes the value of mwni(G). To report the edges involved, we show below how to 
first construct a minimum weight cover of G in 0{^/nW/k(n, W/N)) time and then use 
this cover to construct a maximum weight matching in 0{^/nm/k{n, m)) time. Thus, 
the time required to construct a maximum weight matching is 0{y/nW/k{n, W/N)). 
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Lemma 3.1. Assume that h,GfnCh, and are defined as in §|J. Let he 
any minimum weight cover of G^. If D is a function on V{G) such that for every 
u G V{G), D{u) = Gh{u) + G^{u), then D is a minimum weight cover of G. 

Proof. Consider any edge uv of G. If uv is not in G^, then w{u,v) < Gh{u) + 
Ch{v) < D{u) + D{v). Assume that uv is in G^. Note that its weight in G^ is 
uj{u,v)-Ch{u)-Ch{v). Since G,f is a cover, G^(m)+G^(u) > w{u,v)-Gh{u)-Ch{v). 
Thus, D{u) + D{v) = Ch{u) + C^{u) + Ch{v) + C^{v) > w{u, v). It foUows that D is 
a cover of G. To show that _D is a minimum weight one, we observe that 

E„ey(G)^H - EueviG)CHiu) + C-{u) 

= Y.uev(G)(^hiu) + Y.uev{G)'^hi^) 



mwm(G/i) + mwm(G^) by Fact 2.1 



= mwm(G). by Theorem 2.2 



By Fact p.l| , D is minimum. □ 

By Lemma ^.l| , a minimum weight cover of G can be computed using a recursive 
procedure similar to Compute- MWM as foUows. 
Procedure Compute-Min-Cover(G) 

1. Construct Gi from G. 

2. Find a minimum weight cover Gi of Gi. 

3. Construct Gf from G and Gi. 

4. If Gf is empty, then return Gi; otherwise, let G^ = Compute-Min-Cover(Gf ) 
and return D where for all nodes u in G, = Gi(u) + Gf^(u). 

Theorem 3.2. Compute-Min-Cover(G) correctly computes a minimum weight 
cover of G in 0{y^W/k{n,W/N)) time. 



Proof. The correctness of Compute-Min-Cover(G) follows from Lemma 3.1. For 



the time complexity, the analysis is similar to that of Theorem 2.3. □ 

Now, we show how to recover a maximum weight matching of G from a minimum 
weight cover D of G. 

Procedure Recover-Max-Matching(G, D) 

1. Let H be the subgraph of G that contains all edges uv with w{u, v) = D{u) + 
D{v). 

2. Make two copies of H. Call them iJ" and . For each node u of i?, let 
and u^ denote the corresponding nodes in H"^ and , respectively. 

3. Union and H'' to form and add to the set of edges | u e 
V{H), D{u) = 0}. 

4. Find a maximum cardinality matching K of iJ'^'' and return the matching 
K" = {uv I u'^u'^ e K}. 

Theorem 3.3. Recover-Max-Matching(G, £)) correctly computes a maximum 
weight matching of G in 0{y/nm/k{n,m)) time. 

Proof. The running time of Recover-Max-Matching(G, D) is dominated by the 
construction of K. Since i?"^ has at most 2n nodes and at most 3m edges, K can be 
constructed in 0{^/nm/k{n,m)) time using Feder-Motwani algorithm Q. 

It remains to show that K"^ is a maximum weight matching of G. First, we argue 
that H"-^ has a perfect matching. Let M be a maximum weight matching of G. By 



Fact 2.1, D{u) + D{v) — w{u,v) for every edge uv G M. Therefore, M is al so a 



matching of H. Let U be the set of nodes in H unmatched by M. By Fact |2T 
D{u) = for all u<eU. Let Q be {u"-u'' \ u e U}. Let M" = {u'^w'^ | uu S M} and 
M** = {li^u'' \uv£ M}. Note that Q U Af" U A//'' forms a matching in i/'"'' and every 
node in i?"'' is matched by either Q, A/" or M^. Thus, ff"^ has a perfect matching. 



6 



Kao, Lam, Sung, and Ting 





(a) m 



(b) ^(G)|C, 



Fig. 4.1. (a) The unfolded graph d>( G) of the bipartite graph given in Figure 2.1(a). (b) With 
respect to the cover defined in Figure 2.1(c), the node y\ in <j>(G) is the only node satisfying the 
condition that 1 < Cfi{y). Thus, ij>{G)\Ci^ comprises only the edges incident to yi. 



Since iV' is a maximum cardinality matching of H"''', K must be a perfect match- 
ing. For every node u with D{u) > 0, must be matched by K. Since there is 
no edge between and any in H°'^ , there exists some v°- with u°'v°- G K. Thus, 
every node u with D{u) > must be matched by some edge in K"'. Therefore, 

E™gk° wiu,v) = T,uexuY,D{u)>o^('^) = EuGXur = mwm(G), and K'' is a 

maximum weight matching of G. □ 



4. All-cavity maximum weight matchings. In ^4.1, we introduce the notion 



of an unfolded graph. In §4.2, we use this notion to design an algorithm which, 
given a weighted bipartite graph G and a maximum weight matching of G, computes 
mwm(G — {u}) for all nodes u in G using 0{W) time. 

4.1. Unfolded graphs. The unfolded graph (t>{G) of G is defined as follows. 

• For each node u of G, 0(G) has a copies of m, denoted as , u^, . . . , u", where 
a is the weight of the heaviest edge incident to u. 

• For each edge uv of G, 0(G) has the edges u^v^ ,u'^v^~^ , . . . ,u^v^ , where 
l3 = w(u, v). 



See Figure 4.1(a) for an example. Let M be a matching of G. Consider M as a 
weighted bipartite graph; then, by definition, <j){M) = U™gm{"^^^' ' ' ' i^^w^ | (3 = 
w{u,v)} is a matching of (t>{G). The number of edges in cj){M) is equal to the total 
weight of the edges in M, i.e., |0(M)| = X^tmeA/ '"'(^i ■^)- The next lemma relates G 
and (j){G). 

Lemma 4.1. Assume that M is a maximum weight matching of G. 

1. mwm(G) — mm((/)(G)). 

2. The set 4>{M) is a maximum cardinality matching of (piG). 

Proof. Statement ^ follows from Statement |^. Statement ^ is proved as follows. 
Since M is a maximum weight matching of G, mwm(G) = X^uueM w{u, v) — |(/)(M)| < 
mm(0(G)). By Fact mwm(G) > mm((/)(G)) if and only if mwc(G) > mwc(0(G)). 
We prove the latter as follows. Given a minimum weight cover G of G, we can obtain 
a cover G' of 0(G) as follows. For any node u of G, G'{u^) = 1 if G{u) > and 
i < G(m); otherwise, G'iu^) = 0. Note that w{G') = w(G) = mwc(G). Therefore, 
mwc(G) > mwc(0(G)) and mwm(G) > mm(0(G)). □ 
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4.2. An algorithm for all-cavity maximum weight matchings. Let M be 

a given maximum weight matching of G. 

By Lemma 4.1(2|), (j>{M) is a maximum cardinaUty matching of 4>{G). In light of 
this maximaUty, we say that a path in (/'(G) is alternating for 4>{M) if (1) its edges 
alternate between being in (j){M) and being not in (j){M) and (2) in the case the 
first (respectively, last) node is matched by <j>{M), the path contains the matched 
edge of u as the first (respectively, last) edge. The length of an alternating path is 
its number of edges. An alternating path may have zero length; in this case, the 
path contains exactly one unmatched node. An alternating path P can modify (j){M) 
to another matching, i.e., {(t>{M) UP) — {<j>(M) DP). If P is of even length, the 
resulting matching has the same size as <t>{M). If P is of odd length, P modifies M 
to a strictly smaller or bigger matching; yet the latter is impossible because (f>{M) is 
maximum. Intuitively, we would like to maximize the size of the resultant matching 
and even-length alternating paths are preferred. 

Our new algorithm for computing mwm(G — {u}) is based on the observation 
that mwm(G — {u}) can be determined by detecting the smallest i such that u' has 
an even-length alternating path for (j){M). Details are as follows. 

Definition. For each in (j){G), let p{u^) = if there is an even-length alternating 
path for (j){M) starting from u'; otherwise, let p{u^) = 1. 

The following lemma states a monotone property of p{u^) over different i's. 
Lemma 4.2. Consider any node u in G. Let . .. be its corresponding 

nodes in (f>{G). If p{u^) — 0, then p{u^) — for all j € [i,(3]. Furthermore, there exist 
P ^ i + I node-disjoint even-length alternating paths Pi, Pi+i, ■ ■ ■ Pp for (j){M), where 
each Pj starts fi 



rom u- 



Proof As p(w») = 0, let P, = ul° , vl° ,u1\v{\..., Up'Ti , v^^i , Up" be a shortest 
even-length alternating path for (j){M) where Uq° — u'. 

Based on P^, we can construct an even-length alternating path P^+i for 4>{M) 
starting from as follows. If u^~^^ is not matched by 4>{M), P^+i is simply a 

path of zero length. From now on, we assume that is matched by 0(M). As 
P is of even length, Up'' is not matched by 4>{M). Then, by the definition of 4'{M), 
Up"'^^ is also not matched by 4>{M). Let h be the smallest integer in such 
that u'^^^ is not matched by 0(M). Notice that, for aU i < h, u^^^ is matched 
to w^*~^; furthermore, (/)(G) contains an edge between and Thus, 

Pi+i — w*"*"^, Wq°~^, wj^^"'^, is an even-length alternating path for 

4i{M). Similarly, for j = i-\-2, • • • , /3, we can use Pi to define an even-length alternating 
path Pj for 4){M) starting from By construction. Pi, Pi+i, ■ ■ ■ Pp are node-disjoint. 
□ 

The next lemma is the basis of our cavity matching algorithm. It shows that 
given mwm(G) (i.e., the weight of M), we can compute mwm(G — {w}) from the 
values p(u'); ^-nd the p(m*)'s can be found in 0{W) time. 
Lemma 4.3. 

Si<i</3 = niwm(G) -mwm(G- {u}). 

2. For all € 4>{G), p{u^) can be computed in 0(W) time in total. 
Proof. The two statements are proved as follows. 

Statement |l|. Let k be the largest integer such that p{u'') = 1. By Lemma 4.2 , 
p(m*) — 1 for all 1 < i < fc, and otherwise. Note that if p{u^) — 1, u' must be 
matched by (piM). Thus, X]i<i<;3 /°(^*) = ^- Below, we prove the following two 
equalities: 
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(1) mm(0(G) - {u^,...,u''}) = mm(0(G)) - k. 

(2) mm(0(G) - {u\ . . . , u'^}) = mm(0(G) - {u\ . . . ,u''}). 

Then, by Lemma 4.1, mwm(G) = mm(0(G)) and mwm(G — {u}) = mm(0(G) — 
{u^, . . . , u^}). Thus, mwm(G) — mwm(G — {u}) = k and Statement |l| follows. 

To show Equality (1), let H be the set of edges of 4>{M) incident to with 
1 < i < fc. Let M' = (j>{M)-H. Then, \M'\ = \(t>{M)\-k. We claim that M' is amax- 
imum cardinality matching of (/)(G) — {u^, ...,u'^}. Hence, mwm(0(G) — {u^, u'^}) = 
|(j!)(M)| — k, and Equality (1) follows. We prove the claim by contradiction. Suppose 
M' is not a maximum cardinality matching of (j){G) — {u^, . . . , u*^}. Then, there exists 
an alternating path P that can modify M' to a larger matching of 4>{G) — {u^, . . . , u'^} 
1^,^; in particular, the length of P must be odd and both of its endpoints are not 
matched by M' . P must start from some node with u^v^ € (t>{M) and i < k; 
otherwise, P is alternating for (j){M) in G and (f>{M) cannot be a maximum cardi- 
nality matching of (f>{G). Let Q be a path formed by joining u^v^ with P. Q is an 
even- length alternating path for (t>{M) starting from in <j){G). This contradicts the 
fact that there is no even-length alternating path for 4>{AI) starting from for i < k. 

To show Equality (2), we first note that mm(0(G) — {u^, . . . , u^}) < mm(0(G) — 
{u^ , . . . ,u'^}). It remains to prove the other direction. By Lemma 



4.2 



we can find 

P — k node-disjoint even- length alternating paths Pk+i, ■ ■ ■ jPp for (f>{M), which start 
from • • • , u'^. Pj starts at . Let M" = (0(A/) U (P^+i U • • • U Pp)) - {cj){M) n 

{Pj+i U • • • U Pp)). Note that \M"\ = |0(M)| and there are no edges in M" incident 
to any of • • • , u'^. M" is a matching of 0(G) - {u^+^ , ■■■ ,uf^} and M" - H 

of (j>{G) - {u\ . . .,uf^}. \M" H\ > \M"\ - k ^ \HM)\ - k. Since mm(0(G) - 
{M^ . . .,u''}) = |(/)(M)|-/c by Equality (1), it follows that inm{(f)iG) - {u\ . . .,u''^}) > 
\M" -H\> mm(0(G) - {u\ . . Therefore, Equality (2) holds. 

Statement |[ We want to determine whether p{u^) = for all nodes £ 0(G) in 
0(W) time. By definition, jo(u*) = if and only if there is an even-length alternating 
path for (j){M) starting from m*. Let us partition the nodes of </'(G) into two parts: 
(t){X) = {u^ e 0(G) I w G X} and 0(r) = {u' e 0(G) I u e y}. Below, we give the 
details of computing p{u'^) for all S 0(^). The case where e 0(5^) is symmetric. 

Let D be a directed graph over the node set 4){X). D contains an edge u^v^ if 
there exists a node e 4>(Y) such that u^w^ e 4>{G) — 4>{M) and w^v^ G (l>{M). 
Consider any node of D that is unmatched by 0(Af ). A directed path in D from 
to a node corresponds to a path in 4i{G), which is indeed an even- length alternating 
path for 4){M) starting from w*. Therefore, for any € 4>{X), piu^) = if and only 
if is reachable from some node in D that is unmatched by (j){M). We can identify 
all such u' by using a depth-first search on D starting with all the nodes unmatched 
by M. The time required is 0{\D\). As \D\ < |0(G)| = W, the lemma follows. □ 

The following procedure computes mwm(G — {u}) for all nodes u of G. Let M 
be a maximum weight matching of G. 
Procedure Compute- All-Cavity(G, M) 

1. Construct 0(G) and 0(M). 

2. For every j S [0, n/2], determine Aj from (j){M). 

3. For every node of (p{G), if £ IJ^ then p{u^) ~ 0; otherwise p{u^) = 1. 

4. For every node u of G, compute mwm(G — {u}) = mwm(G) — J2i<i<p 
where u^, v?, . . ., are the nodes corresponding to u in (p{G). 

Theorem 4.4. Compute- All-Cavity(G, M) correctly computes mwm(G — {u}) 
for all u of G in 0{W) time. 

Proof. Follows from Lemma 4.3 □ 
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